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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 
RESEARCH  MEMORANDUM 

THEORETICAL  STUDY  OF  THE  TUNNEL-BOUNDARY  LIFT 
INTERFERENCE  DUE  TO  SLOTTED  WALI£  IN  THE 
PRESENCE  OF  THE  TRAILING -VORTEX 
SYSTEM  OF  A LIFTING  MODEL 
By  Clarence  W.  Matthews 

SUMMARY 


The  equations  presented  in  this  paper  give  the  interference  on  the 
trailing-vortex  system  of  a uniformly  loaded  finite- span  wing  in  a 
circular  tunnel  containing  partly  open  and  partly  closed  walls,  with 
special  reference  to  symmetrical  arrangements  of  the  open  and  closed 
portions.  Methods  are  given  for  extending  the  equations  to  include 
tunnel  shapes  other  than  circular.  The  rectangular  tunnel  is  used  to 
demonstrate  these  methods.  The  equations  are  also  extended  to  non- 
uniformly  loaded  wings. 

An  analysis  of  the  equations  for  certain  configurations  has  shown 
that:  (l)  only  a small  percentage  of  slot  opening  is  required  to  give 

zero  interference  conditions  if  the  tunnel  contains  four  or  more  slots; 
(2)  in  the  configurations  studied,  the  ratio  between  the  slotted-tunnel 
interference  and  the  closed- tunnel  interference  at  the  center  of  the 
tunnel  is  approximately  constant  for  various  model  spans;  and  (3)  tunnels 
containing  an  odd  number  of  slots  or  nonsymmetrical  slot  arrangements 
cause  an  additional  rolling  moment  or  a cross  flow  on  the  wings,  or  both. 


INTRODUCTION 


In  a study  of  solid -blockage  interference  (see  ref.  l),  it  has  been 
shown  that  tunnels  containing  mixed  boundaries,  that  is,  partly  open  and 
partly  closed  walls,  will  eliminate  or  greatly  reduce  such  interference. 
Such  tunnels  have  been  shown  to  be  very  useful  test  equipment  (ref.  2). 
Since  the  slotted  tunnel  configuration  required  to  eliminate  solid 
blockage  may  not  eliminate  lift  interference,  it  is  necessary  to  study 
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Since  the  slotted  tunnel  configuration  required  to  eliminate  solid 
blockage  may  not  eliminate  lift  interference,  it  is  necessary  to  study 
the  interference  on  the  trailing  vortices  of  a lifting  model  in  order 
to  make  the  necessary  corrections  to  the  lift  characteristics  of  a model. 

The  problem  of  one  or  two  slots  has  been  treated  by  various  authors 
(see  refs.  3 to  6).  The  case  of  more  than  two  slots  has  also  been  treated 
in  references  7 and  8.  Reference  7 treats  only  small  wings  in  circular 
tunnels,  and  reference  8 treats  only  the  case  for  large  numbers  of  evenly 
spaced  slots. 

The  purpose  of  this  report  is  to  present  equations  which  express 
the  tunnel-wall  interference  due  to  mixed  open  and  closed  boundaries  in 
the  presence  of  the  trailing- vortex  system  of  a finite-span  lifting  model 
at  subsonic  velocities.  Special  attention  is  given  to  test  sections  in 
which  the  slots  are  symmetrically  located  with  respect  to  both  axes. 

Various  extensions  of  the  theory  have  been  made,  following  in  a 
general  fashion  the  methods  of  reference  7.  These  extensions  include 
the  effects  of  wing  span,  slot  configurations,  interference  at  points 
near  the  center,  nonuniform  loading,  and  methods  for  calculating  the 
interference  in  tunnels  of  other  than  circular  cross  section. 

Numerical  calculations  of  the  interference  characteristics  of  several 
symmetrical  cases  are  presented  and  are  used  to  show  the  properties  of 
the  interference  of  circular  tunnels  containing  1,  2,  4,  8,  and  12  slots 
symmetrically  located  with  respect  to  the  x-  and  y-axes,  and  a square 
tunnel  containing  8 slots  symmetrically  located  in  the  top  and  bottom 
walls. 


An>^n 

an>bn 
an>  3n 
b 
C 

CL 


SYMBOLS 

> constants 

half -span  of  wing 

tunnel  cross-sectional  area 

lift  coefficient  of  model 


CONFIDENTIAL 


SECURITY  INFORMATION 


NACA  RM  L53A26 


CONFIDENTIAL 


3 


ACd 

ds- 

g 

k 


1 

m 

M 


n,s,q 


P 
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R.P. 


I.P. 


S 

u 

v 

V 


z 

dw 

dz 

r 

7 


drag  of  model 

drag  increment  due  to  interference 
tangential  line  element 
multiple-product  index 

quality  factor,  the  ratio  at  any  given  point  of  the  inter- 
ference of  a slotted  tunnel  to  the  interference  of  a 
closed  tunnel  with  the  same  cross  section 

half -span  length  of  wing  in  £ -plane 

number  of  slots  in  tunnel 

Mach  number 

summation  indices 

special  number  defined  in  equation  (2) 

radial  distance  of  point  from  coordinate  center 

real  part  of  complex  function 

imaginary  part  of  complex  function 

model  or  wing  area 

x-component  of  velocity 

y-component  of  velocity 

velocity  vector 

coordinates  of  rectangular  system  (fig.  1) 
complex  variable,  x + iy 

complex  velocity,  u - iv 

circulation  about  a point,  positive  in  counterclockwise 
direction 

function  of  0n  and  b defined  in  equation  (32) 
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a 

t 

e 

$ 

\ 

<P 


proportionality  factor  » (see  ref.  9) 

interference  correction  angle  to  be  added  to  measured 
angle  of  attack 

strength  of  a source 

complex  variable  in  transformed  plane 

angular  coordinate  of  polar  system  (fig.  1) 

complex  function  of  z (z  = el13) 

parameter  defined  in  equation  (105)  (also  see  par.  14.8, 
ref.  10) 

perturbation  potential 


dcp 

8n 


derivative  of  potential  in  direction  normal  to  a given  line 


THEORY  OF  LIFT-VORTEX  INTERFERENCE  IN  SLOTTED  TUNNELS 

General  Analysis 


Theoretical  boundary  conditions  of  flov  about  trailing  vortices  in 
a tunnel  containing  mixed  open  and  closed  portions.-  The  equations  for 
the  interference  on  the  lift  of  a model  due  to  mixed  open  and  closed 
tunnel  walls  in  the  presence  of  a trailing- vortex  system  may  be  obtained 
by  considering  the  same  two-dimensional  approximation  of  the  flow  field 
that  is  used  in  reference  9»  The  conditions  of  this  two-dimensional 
approximation  may  be  briefly  stated  as:  (l)  the  tunnel  and  its  boundaries 

extend  from  a point  an  infinite  distance  upstream  of  the  model  to  a point 
an  infinite  distance  downstream  of  the  model;  (2)  the  velocities  induced 
by  the  trailing  vortices  in  the.  cross-sectional  plane  located  at  the  model 
are  one-half  of  those  induced  in  a far-downstream  cross-sectional  plane; 
(3)  the  induced-velocity  flow  field  in  the  far-downstream  section  may  be 
treated  with  two-dimensional  methods;  (4)  the  boundary  condition  which 
must  be  satisfied  at  a solid  portion  of  a tunnel  wall  is  that  the  flow 
must  be  tangential  to  it,  or 


I®  = 0 (1) 

on 
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(5)  the  condition  which  must  be  satisfied  at  an  open  portion  of  a tunnel 
wall  is  that  the  potential  over  that  portion  must  be  constant  or  the  flow 
must  be  normal  to  it;  (6)  no  singularities  other  than  the  trailing  vortices 
can  exist  within  the  boundaries  of  the  tunnel;  and  (7)  the  constant 
potential  in  every  slot  must  be  equal  to  that  in  every  other  slot.  The 
final  condition  is  required  because  the  pressure  is  the  same  at  every 
slot  and  may  be  shown  by  the  following  considerations.  Since  the  pres- 
sure is  constant  over  the  entire  region  outside  the  tunnel  there  will 
be  no  pressure  differences  between  the  slots  at  a point  far  upstream, 
and  hence  there  can  be  no  flow  between  the  slots  due  to  external  influ- 
ences. Also,  since  this  point  is  too  far  from  the  model  to  be  influ- 
enced by  it,  there  can  be  no  flow  due  to  the  model.  As  there  is  no  flow 
between  the  slots,  no  potential  gradient  can  exist  between  the  slots; 
therefore,  the  potentials  in  all  the  slots  must  be  the  same.  It  there- 
fore follows  by  condition  (5)  that  in  the  far-downstream  position,  the 
constant  potential  in  every  slot  is  equal  to  that  in  every  other  slot. 

Velocity  fields  in  circular  slotted  tunnels.-  The  previously  stated 
boundary  conditions  may  be  satisfied  by  using  complex  velocity  functions 
rather  than  complex  potential  functions.  This  is  done  by  selecting  a 
complex  velocity  function  that  has  singularities  at  the  wing  tips  and 
has  a flow  direction  either  normal  or  tangential  to  the  tunnel  walls. 

If  this  function  is  multiplied  by  another  complex  function  whose  value 
on  the  tunnel  wall  changes  from  all  real  to  all  imaginary,  or  the 
opposite,  at  each  slot  edge,  then  the  flow  will  be  rotated  90°  at  those 
points  so  that  the  final  flow  of  the  product  of  the  two  functions  will 
be  normal  to  the  wall  on  selected  portions  and  will  be  tangential  to 
the  wall  on  the  remaining  portions.  It  may,  however,  be  expected  that 
the  second  function  will  introduce  singularities  within  the  tunnel  which 
are  not  permitted  according  to  the  stated  boundary  conditions,  so  that  a 
third  function  which  contains  all  the  forbidden  singularities  must  be 
used  in  such  a fashion  that  it  will  cancel  the  forbidden  singularities 
of  the  second  function. 

In  order  to  make  up  the  first  function,  suppose  that  the  two  singu- 
larities at  the  wing  tips  with  their  reciprocal  singularities  are  written 

as  l/(z2  - b2)(l  - z2b2),  where  z is  the  complex  variable  x + ly 
and  b is  the  semispan  of  the  lifting  wing.  If  this  function  is  examined 
by  letting  z = e1-®,  where  0 is  a polar  angle  (fig.  1),  it  will  be  found 
that  the  flow  may  be  made  normal  to  the  walls  if  the  factor  z is 
included  in  the  numerator.  This  flow  may  also  be  made  tangential  to  the 
wall  by  multiplying  by  the  factor  i.  Thus,  the  first  function  may  be 
written 


pz 

(z^  - b2)^l  - z2b2) 


(2) 
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The  symbol  p may  be  chosen  to  be  either  1 or  i,  depending  on  whether 
normal  or  tangential  flow  is  required  at  the  walls. 


The  second  function  may  be  developed  by  considering  the  square  root 
of  a function  which  is  real  on  the  wall  and  changes  sign  at  each  slot  edge 
so  that  the  square  root  of  the  function  changes  from  real  to  Imaginary  at 


each  slot  edge, 
where  z = e ^ 


Such  a function  may  be  expressed  by  \/ cot  — - cot  ^ 

2 2 

and  i3  is  in  general  complex.  Examination  of  this  func- 


{■ 


e 

tion  shows  that  it  becomes  wcot  ^ - cot  on  the  tunnel  wall,  where 
0g  is  the  polar  angle  of  a slot  edge  (see  fig.  1).  The  term  under  the 


radical  changes  from  positive  to  negative  at  0 = 0g  so  that  the  func- 
tion changes  from  real  to  imaginary  at  the  point  0g  on  the  tunnel  wall. 

If  a number  of  these  functions  having  different  0g's  marking  the  transi- 
tion from  open  to  closed  sections  of  the  tunnel  wall  are  multiplied 
together,  it  may  be  seen  that  the  product  will  change  sign  at  each  value 
of  0rr  so  that  the  function  will  be  real  on  alternate  sections.  This 

O 

use  of  cot  & also  suggests  the  use  of  other  complex  trigonometric  func- 

u 

tions  such  as  cos  fl,  sin  fl,  and  tan  ^ . The  sine  and  cosine  terms  will, 

on  examination,  be  found  to  introduce  two  slot  edges  for  each  value  of  0g, 
rather  than  one.  These  two  slot  edges  will  be  found  to  be  located  symmet- 
rically with  respect  to  the  x-axis  for  the  cosine  term  and  to  the  y-axis 
for  the  sine  term.  Thus,  the  use  of  these  functions  is  suggested  for 
symmetrical  slot  configurations. 


Since  any  of  these  functions  may  introduce  singularities  within  the 
tunnel,  it  is  necessary  to  examine  them  for  such  singularities.  The  four 
functions  may  be  written  as  the  following  multiple  products. 


n icos  $ - COS  0g 
g=l 


^COS  d - COS  ©1  ^COS  - cos  0g 


/COS  d - COS  0 


6 


■d  - sin  0g 


(3) 

(4) 


<1 


n 

g-i 


icot  — - cot  — 
12  2 


(5) 
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and 

I!  ^tan  | - tan  -1  (6) 

When  z = eil3  is  substituted  for  $,  the  above  equations  become, 
respectively, 


2-q/2z-q/2 


(7) 


2-q/2z-q/21-q/2  g p . I - 2iz  sin  0g  (8) 


(z  - l)'^2  n /i(z  + 1)  - (z  - l)cot  ^ (9) 

g=l 


(z  + u-a-*/*  n i/(z  - 1)  - + i)taa  5 d°) 

g=l  v 2 

Examination  of  the  functions  (7)  to  (10)  shows  that  they  contain 
forbidden  singularities  at  the  following  values  of  z:  the  cosine  and 

sine  terms  at  z = 0,  the  cotangent  term  at  z = 1,  and  the  tangent 
term  at  z = -1.  Several  functions  which  contain  singularities  identical 
to  those  appearing  in  functions  (7)  to  (10)  and  which  are  real  on  the 
tunnel  are 


(an  cos  ntf  + bn  sin  n$) 


(11) 
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which  has  a singularity  at  z ■ 0} 

(on  + i&^cot^l)  (12) 

which  has  a singularity  at  z = lj  and 

+ i3n)tan^j  (13) 


which  has  a singularity  at  z = -1;  where  an,  an,  bn,  and  3n  are 
real  constants.  It  may  be  seen  now  that  by  dividing  the  trigonometric 
series  by  the  multiple  product  which  contains  the  same  singularities, 
the  forbidden  singularities  will  be  canceled  out  of  the  final  equation, 
leaving  an  equation  which  has  only  the  desired  singularities  within  the 
tunnel. 


The  reason  that  the  multiple  product  is  chosen  for  the  denominator 
may  be  seen  by  examining  the  flow  about  the  slot  edges.  Since  the  flow 
must  turn  a sharp  corner  as  it  goes  around  the  edge  of  each  slot,  the 
velocity  must  be  infinite  at  that  point.  Since  the  multiple-product 
function  becomes  equal  to  zero  at  each  slot  edge,  it  must  be  placed  In 
the  denominator  of  the  final  complex  velocity  function  to  insure  the 
required  infinite  velocities. 


Before  writing  the  final  complex  velocity  function,  it  is  observed 
that  function  (2)  has  a zero  at  the  point  z = 0 which  must  also  be 
removed  by  either  including  an  extra  term  in  summation  (ll)  or  ty 
multiplying  summation  (12)  or  (13)  by  a^  cos  3 + b^_  sin  -3.  With  this 
note  in  mind,  the  final  complex  velocity  functions  may  be  written: 


dw 

dz 


pz  (an  cos  ntf  + bn  sin  ntf) 

n=0 


(z2  - b2)(l 


(U) 
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dw 

3z  = 


(an  cos  nd  + bn  sin  nd) 


q. 

- z^b2 ) J J ^fsin  ^ - sin  0g 
g=l 


(15) 


dw 

3?  = 


pz(a]_  cos  -a  + bi  sin  d) 


q/2  . v 

rj  (an  + i3n)cot°^-J 


- cot  !fi 
2 


(16) 


and  an  equation  similar  to  (l6),  using  tan  ^ rather  than  cot 

In  the  application  of  these  equations,  it  may  be  observed  that  each 
value  of  0g  introduces  two  slot  edges  into  each  of  equations  (l4) 
and  (15)  and  only  a single  slot  edge  into  equation  (l6).  Since  each 
panel  has  two  edges,  q must  be  equal  to  the  number  of  panels  m in 
equations  (lk)  and  (15)  and  to  twice  the  number  of  panels,  or  2m,  in 
equation  (l6).  Also,  if  equation  (14)  or  (15)  Is  applied  to  tunnels 
containing  an  odd  number  of  slots,  it  will  be  found  that  the  singularity 
arising  from  the  multiple  product  will  contain  a term  of  the  order  l/2, 
which  cannot  be  removed  by  the  summation.  In  order  to  remove  this  singu- 
larity, the  summation  of  equation  (l4)  or  (15)  must  be  rewritten  as 

m+1 

5^  (*n  cos  2-n-2+-  -1  d + bn  sin  2n- £ 1 d)  (17) 

This  series  will,  upon  examination,  be  found  to  contain  a singularity  of 
the  order  l/2,  which  will  cancel  the  one -half -power  term  due  to  the  exten- 
sion of  the  multiple  product  over  an  odd  number  of  slots. 

In  order  to  make  the  final  application  of  equations  (l4),  (15), 
or  (16)  to  a wind  tunnel,  the  arbitrary  constants  of  the  suusnations  must 
be  evaluated.  This  evaluation  can  be  effected  by  using  the  remaining 
boundary  conditions,  which  are  that  the  potentials  in  each  slot  must  be 
equal  to  each  other  and  that  only  a vortex  flow  can  exist  about  the 
singularities  within  the  tunnel. 


CONFIDENTIAL 


10 


CONFIDENTIAL 


NACA  RM  L53A26 


The  potential  condition  may  "be  evaluated  by  the  following  line 
integral: 


P(x 2J2)  _ 

/ V*ds  = cp/  s - cp,  , = 0 (l8) 

J(X;  i,yi)  (X2’y2) 

in  which  the  limits  of  integration  terminate  in  the  slot.  One  path  which 
may  be  used  is  the  streamline  which  flows  along  the  panel.  'When  this 
streamline  is  used  as  the  path  of  integration,  equation  (l8)  becomes 


re 


n+1 


v d0  = <Pn+l  - «Pn  = 0 


(19) 


as  the  potential  is  the  same  over  all  slots.  Since  the  velocity  V will 
in  general  be  complex,  two  equations  exist  for  each  panel. 

The  strengths  of  the  circulation  about  a pole  and  of  the  source  at 
a pole  are  fixed  by  the  relation, 


r + io  = £ — dz  (20) 


where  r is  the  circulation  due  to  lift  and  a is  the  source  strength 
and  is  equal  to  zero.  Since  equation  (20)  must  be  evaluated  about  each 
pole,  it  gives  four  equations  which  may  be  used  to  evaluate  the  constants. 

The  set  of  boundary  equations  (19)  and  (20)  may  be  considered  as 
a set  of  simultaneous  equations  in  the  unknown  constants  found  in  equa- 
tion (14),  (15),  or  (16)  and  is  used  to  determine  the  values  of  these 
constants.  If  these  constants  are  to  be  uniquely  determined,  there  must 
be  as  many  equations  as  there  are  constants.  An  examination  of  equa- 
tions (19)  and  (20)  shows  that  equation  (19)  gives  2m  equations  and 
that  equation  (20)  gives  four  equations,  making  a total  of  2m  + 4 
equations  which  may  be  used  to  determine  the  unknown  constants.  The 
number  of  constants  which  must  be  satisfied  is  determined  frcm  an  exami- 
nation of  equation  (l6).  It  is  shown  in  the  paragraph  following  equa- 
tion (l6)  that  q must  equal  2m,  so  that  the  number  of  constants  must 
equal  2m  + 4;  hence  the  number  of  equations  and  the  number  of  constants 
are  equal,  so  that  all  the  constants  may  be  uniquely  determined,  thus 
giving  the  problem  a unique  solution. 
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Derivation  of  Equations 


Slots  symmetrically  located  with  respect  to  the  x-  and  y-axes.- 
In  the  case  of  symmetry  about  the  x-  and  y-axes,  considerable  simplifica- 
tion results.  The  simpler  equation  (14)  may  be  used  in  place  of  equa- 
tion (l 6).  In  equation  (14),  it  may  be  observed  that  each  value  of  0g 
in  the  multiple  products  can  be  used  to  produce  four  changes  of  sign, 
provided  these  products  are  written  as 

m/2 

n/  COS2iS  - COS^0g  (21) 

6=1 


m/2 

n /sin2 3 - sin20g  (22) 

g=l 


These  changes  of  sign  are  observed  at  0g,  -0g,  * - 0g,  and  jt  + 0g, 
or  one  in  each  quadrant.  Hence,  if  the  values  of  0g  are  known  in  one 
of  the  quadrants,  the  multiple  product  need  be  extended  only  over  that 
quadrant,  as  all  the  other  slots  will  be  automatically  introduced.  If 
the  slots  are  evenly  spaced  and  of  equal  width,  function  (21)  or  (22) 
can  be  replaced  by 


or 


| * 

\p  \ * 


- COS2  — 0t 
2 1 


- sin2  S 0 


2 


(23a) 


(23b) 


That  (23a)  or  (23b)  produces  slots  in  the  desired  location  may  be  seen 
by  substituting 


- Hus  + 

m ^ 


on  the  tunnel  wall  and 


0-,  = — - 
x m 


qPl 


where  m is  the  number  of  slots  in  the  tunnel  and  n is  any  integer. 
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Taking  this  modification  into  consideration,  equation  (14)  may  be 
written  with  the  substitution  of  z = e^3  as 


5 + 1 


pz 


dw 

dz 


*n 


s2n  + 1 
2zn 


- i*n 


s2n  - lN 
2zn  ; 


m, 


■/2 


(24) 


z“  - 


b2)  (l  - A2)-  - 1 , n 1 - 2z2cos  26.  + 1 
2m/2zm/2  g=l  ' S 


which  may  be  reduced  to 


m , 

l-1  P 

22  iT 


dw 

dz 


m . 
7T-  n+1 


m 


n+1 


anz2  (z2n  + l)  - ibnz2  (z2n  - l) 

n=0  L_ 


m/2 


(25) 


(z2  - b2)  (l  - z2b2)  "[  [ ^ - 2z2cos  20g  + 1 

g=l 


where  the  factor  p is  equal  to  1 if  slots  intersect  the  x-axis  and  is 
equal  to  i if  panels  intersect  the  x-axis.  If  the  slots  are  evenly 
spaced,  and  of  equal  width,  the  function 


Jz2*  - 2zmcos  mQ^.  + 1 

may  be  substituted  for  the  multiple  product  in  equation  (24)  or  (25). 

The  constants  of  equation  (25)  must  be  evaluated  by  using  the  bound- 
ary condition  for  the  equality  of  the  potential  in  each  slot,  equa- 
tion (19) , and  for  the  circulation  and  nonexistence  of  sources,  equa- 
tion (20).  Equation  (20)  is  evaluated  by  applying  Cauchy's  integral 
theorem  to  equation  (25)  to  obtain  the  required  line  integral  of  dv/dz, 
or 


T + ic  = 


m 1 _ 

m . o + 1 r E i 

2nip22  5Z  ki*2  n+  (*>2n  + l) 
n=0  I— 


m - 

ibnb2  ' n+  (b2n  - l) 


m/2 


(26) 


2b  (l  - b4)  J1  /b4  - 2b2cos  20g  + 1 
g=l 
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In  solving  equations  (26)  to  (28)  for  an  and  bn,  it  is  observed 
that  for  each  a^: 


an  ~ ^ 


y.l1  -.^1  n fk  - 2b2 


*2 


|-1 


cos  20g  + 1 


(29) 


with  a corresponding  relation  between  bn  and  a new  constant  Pn.  With 
the  substitution  of  these  new  constants  On  and  pn,  equation  (25)  may 


be  written 


m/2 


dw 

dz 


Pbr(i  - b4)  n - 2b2cos  20^.  + 1 ^+1 

g=l 


7T 


jt(z2  - b2)  (l  - z2b2)  T~[  \jz^  - 2z2cos  2 + 1 n 

g=l 


m 

n+1 


(z2n  + l) 


=0 


5-n+l 


iPn2' 


(,2n  - 1) 


and  equations  (26),  (27),  and  (28)  become,  respectively, 


m 


+ 1 


i 

P 


£-n+l 


— - n+1 


(b2n  + l)  - IB,*2  (b2n  - l) 


( j0) 


( jla) 


" p = ll  Xn(-b)?”n+1l(-b)2n  + l]  - iPn(-b)2  "n+1[(-b)2n  - l]| 


(31b) 
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■/ 


where  <%  and  3n  are  the  new  constants. 

It  can  he  seen  by  examining  equations  (31)  that  a number  of  the 
a^'s  and  3n's  are  equal  to  zero.  This  is  shown  by  first  considering 
the  nature  of  the  simultaneous  equations  (31).  The  equating  of  the  real 
and  imaginary  parts  divides  the  entire  set  into  two  separate  sets  of 
simultaneous  equations,  one  of  which  involves  a^'s  only  and  the  other 
3n's  only.  When  p is  chosen,  one  of  the  sets  will  be  homogeneous 
whereas  the  other  set  will  contain  two  equations  with  constant  terms 
equal  to  -1.  Thus,  one  set  of  constants  becomes  equal  to  zero  whereas 
the  other  set  will  have  values  which  may  be  different  from  zero.  This 
remaining  set  of  equations  can  be  reduced  by  considering  the  values  of 
the  integrals  across  any  one  of  the  panels  and  its  counterpart  which  is 
symmetrically  located  with  respect  to  the  y-axis.  It  can  be  shown  by 
substituting  it  - cp  for  0 in  these  integrals  that  they  will  have  the 
same  absolute  values  and  the  same  sign  if  n is  even  but  opposite  signs 
if  n is  odd.  The  same  feature  is  also  observed  in  comparing  equa- 

tlons  (31a)  and  (31b):  that  ls-  tf  I - n + 1 ls  odd’  then  the  tenns 

which  contain  b to  that  power  will  have  opposite  signs  whereas  the 
remaining  terms  will  have  the  same  signs.  It  can  now  be  seen  that,  by 
adding  or  subtracting  equations  (31a)  and  (31b)  and  the  equations  for 
the  integrals  across  each  of  the  symmetrically  located  panels,  it  is 
possible  to  reduce  the  set  of  equations  in  (or  in  3n)  to  two  new 
sets,  one  of  which  contains  the  odd  values  of  n and  the  other  the  even 
values  of  n.  Again  one  of  these  sets  is  homogeneous,  and  one  equation 
of  the  other  set  has  a constant  term  equal  to  -1  so  that  all  the  constants 
associated  with  either  the  odd  values  of  n or  the  even  values  of  n 
will  be  equal  to  zero  and  the  other  set  will  have  values  which  can  be 
different  from  zero. 

This  analysis  shows  that  one  of  four  sets  of  constants  <xgs,  ag8+1, 
3gs»  and  32s+l  can  occur>  thus  suggesting  four  different  synmetrical 
slot  configurations.  The  possible  synnetrical  configurations  which  can 
occur  are  (I)  panels  intersecting  both  axes,  (II)  slots  intersecting  both 
axes,  (III)  a panel  intersecting  the  x-axis  and  a slot  intersecting  the 
y-axis,  and  (IV)  a slot  intersecting  the  x-axis  and  a panel  intersecting 
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m , 
2 + 1 


^ (an  cos  n0  - i3n  sin  n0)d0 
n=0 


(1- 


m/2 


(31c) 


2b^cos  20  + b^)  _[_  J_  /cos^0  - cos^0| 
g=l 


g 
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the  y-axis.  These  conditions  are  found  to  occur  when  the  tunnel  has 
4s  slots  with  p equal  to  i (case  I)  or  1 (case  II)  or  2 (2s  + 1) 
slots  with  p equal  to  i (case  III)  or  1 (case  IV).  Since  there  are 
four  sets  of  constants  and  four  different  symmetrical  slot  configura- 
tions, it  is  to  he  expected  that  each  set  of  constants  can  he  associated 
with  one  of  the  symmetrical  slot  configurations.  When  equations  (31) 
are  set  up  for  each  of  the  symmetrical  slot  configurations,  it  is  found 
that  the  following  associations  exist. 

Case  I,  panels  intersecting  hoth  axes:  the  °l(2s+1)  set  has  va^ues 

different  from  zero. 

Case  II,  slots  intersecting  hoth  axes:  the  P(2s+1)  8et  ^as  val-ues 

different  from  zero. 

Case  III,  panel  intersecting  the  x-axis  and  slot  intersecting  the 
y-axis:  the  set  has  values  different  from  zero. 

Case  IV,  slot  intersecting  the  x-axis  and  panel  intersecting  the 
y-axis:  the  P2s  se"t  has  values  different  from  zero. 

It  is  now  possible  to  write  equations  (30)  and  (31)  by  using  only 
the  constants  which  may  have  values  other  than  zero.  However,  equa- 
tions (30)  and  (31)  can  he  written  in  a more  symmetrical  fonn  if  the 
solution  is  carried  a step  further.  The  equations  for  the  set  of  con- 
stants that  are  used  for  a given  configuration  may  he  written,  using 
the  set  a^s+i  as  an  example, 


-i  = + a3703  + 

0 - Vll  + a3713  + • 
0 = 0^721  + a3?23  + * 


’ “2s+1^0,2b+1 
“28+171,28+1 


“23+172,28+1 


> 


(32) 


where  the  7's  are  the  corresponding  functions  of  9n  and  h.  Now 
consider  as  a matrix  the  coefficients  of  the  right-hand  side  of  the 
equations  whose  constant  terms  are  equal  to  zero,  and  let  A^,  A^, 

A5,  ...  he  equal  to  the  determinants,  respectively,  which  remain 
after  the  column  which  corresponds  to  the  number  2s  + 1 of  the  con- 
stant A28+i  iB  removed.  Then  using  Cramer's  rule, 
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al 


^1 

Al?0±  ~ a3^03  + • • • (-1)Sa2s+1^o,2s+1 


> 


(33) 


ou  = 


Voi  “ V03 


(-1)Sa2s+1^0,2s+1 


with  corresponding  equations  for  the  remaining  values  of  ags+1 

When  the  solution  for  each  constant  equation  (33)  is  substituted 
into  equation  (30),  it  may  be  written 


dw 

37  ~ 


m/2  

-pbr(l  - bM  1 [ yb^  - 2b2cos  2 9 + 1 

g=l  S 

m/2 

it(z2  - b2)  (l  - z2b2)  j_  J_  - 2z2cos  2 0g  + 1 

g=l 


(3D 


The  above  equation  must  be  multiplied  by  one  of  the  following  terms,  the 
choice  depending  upon  the  symmetry  of  the  desired  slot  configuration; 
for  case  I,  or  panels  intersecting  both  axes, 


f (-DSa2s+1z2"ZSE2(2S+1)  + l] 

15=0 

f (-DsA23ti^‘2sE2(2s+1)  ♦ 1] 

s=0 


(35a) 


for  case  II,  or  slots  intersecting  both  axes  (note  that  in  this  case  as 
well  as  in  case  IV,  B has  the  same  relation  to  0 as  A has  to  a 
in  the  demonstrated  case), 

if  (-dws‘2sE2(2s+i)-E1 

fa=1J  — — — (35b) 

28  |^2(2s+1)  _ 1 


'k  5 

("i)S®2s+l^2 


CONFIDENTIAL 


18 


CONFIDENTIAL 


NACA  RM  L53A26 


for  case  III,  or  a panel  intersecting  the  x-axis  and  a slot  intersecting 
the  y-axis, 


m+2 


m+2 

T 


<-i)%J-2stl(>  * i) 


m 


- 2s+l 


s=0 


^ (-D^b2  (b4S  + J 


(o5c) 


and  for  case  IV,  a slot  intersecting  the  x-axis  and  a panel  intersecting 
the  y-axis, 


m+2 


s=0 


DSB 2s^ 


m O 1 

o - 2s+ 1 


(■>  - l) 


m+2 

? (-+)SB2Sb: 

s=0 


| -2S+1,  . 

2 L 4s 


- 1 


(j,5d) 


The  constants  A2S+i,  B2s+1’  Ags,  811(1  ®2s  are  determined  from  the 
following  matrices  in  the  manner  discussed  in  the  material  following 
equation  (32).  The  matrix  for  case  I is 


r,0l 

cos  0 dO 

"°i 

-°1 

fq(0) 

J“°l 

r03 

cos  0 d0 

°2 

fx(o) 

cos  3 0 dO 

o)"'" 


'0m 

2'1 


J°»  « 
2"  c 


cos  0 do  dO 


fj.(0) 


cos^S  + dO 

fl(0) 


(j5e) 
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where 


m/2  

f1(9)  « (l  - 2h2cos  20  + h1*)  IT  )/cos20  - cos20g 


The  matrix  for  case  II  is 


■^0, 

.01 


2 sin  0 d0 

~1^qT 


J0 


J03 


sin  0 d0 

T2 W 


'm/2 


l-1 


sin  0 d0 

tz^) 


02  sin  30  d0 

i f2U) 


(35f) 


where 


m/2  _ 

f (g)  = (l  - 2h2cos  20  + "b  V T"T  Vc°s20g  " cos  9 
2 g=l 
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The  matrix  for  case  III  is 


cos  2 0 d0 
fl(0) 


(35g) 


where 


m/2 

fx(0)  = (l  - Zb2cos  20  + b 7 fj  /cos20  - cos20g 

g=l 


The  matrix  for  case  IV  is 

do  , f0  2 

sin  20  d.0 

jg  f2(e)  J9l 


W sin  20  d0 

f2^0) 


J03 


>UI  <| 
2" 


l-2 


sin  40  d0 

f2(0) 


sin  20  d.0 

f*2(0) 


J6! 


,02  sin(^  + lj  0 d0 

f2(0) 


(35h) 
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where 


m/2 

£2(0)  = (l  - 2b2cos  20  + b^)  I _[  j^cos^e^- -~cos^0 

g=l 


Solution  for  symmetrically  loaded  wings.-  Syranetrically  loaded  wings 
can  be  assumed  to  be  made  up  of  vortex  pairs  of  circulation  -dr  and  dr. 
Since  the  circulation  can  be  expressed  as 


r = r0f (x)  (36) 

the  strength  of  each  pair  becomes 

dr  = rof  *(x)dx  (37) 

This  value  of  the  circulation  may  be  substituted  into  equation  (34)  to 
give  the  contribution  of  each  elemental  vortex  to  the  total  flow.  It  is 
then  necessary  to  integrate  the  equation  over  the  entire  span  to  obtain 
the  flow,  or 


where  (dw/dz)n  is  the  complex  velocity  of  the  nonuniformly  loaded  wing, 
and  (dw/dz)u  is  the  complex  velocity  (eq.  (3*0)  of  the  uniformly  loaded 
wing  with  x substituted  for  b and  a circulation  of  unity. 

Corrections  for  Interference  due  to  lift.-  The  interference  complex- 
velocity  field  at  the  far-downstream  position  is  determined  by  subtracting 
the  complex  velocity  of  the  free  field  from  the  complex  velocity  of  the 
constrained  field,  or 


dwi  _ dw  ibT 

dz  ~ dz  + «(z2  - b2) 


(39) 


A useful  parameter  which  indicates  the  effect  of  slotting  a tunnel 
is  the  ratio  of  the  interference  of  that  tunnel  to  the  interference  of 
a closed  tunnel.  This  parameter,  which  is  the  same  for  the  far-downstream 
position  as  at  the  model  and  which  will  be  called  the  quality  factor  k, 
may  be  expressed  mathematically  as  follows: 
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k = 


dw  ibf 

dz  * nlz2  - t2) 


dwc  ib  r 

dz  + *(Z2  - b2) 


(40) 


where  dwc/dz  is  the  complex  velocity  of  the  closed  tunnel  configuration. 


Once  the  values  of  k are  determined  as  a function  of  the  semi span  b 
for  a specific  tunnel,  the  interference  may  be  calculated.  It  is  shown 
in  reference  9 that  the  interference  on  the  lift  may  be  expressed  as  an 
angle  which  may  be  added  to  the  measured  angle  of  attack  to  obtain  the 
true  or  free-f light  angle  of  attack.  Reference  9 gives  the  correction 
angle  in  radians  as 


eI  = 


&scL 

c 


(41) 


where  5 is  a number  which  is  determined  from  the  geometry  of  the  tunnel 
configuration  and  S is  the  wing  area,  C the  tunnel  area,  and  CL  the 
lift  coefficient.  Also  from  reference  9,  the  increment  which  must  be 
added  to  the  measured  drag  to  obtain  the  correct  drag  is 

qp  2 

acd  = 5 -jr-  (42) 


In  the  case  of  the  slotted  tunnel,  it  is  convenient  to  use  the  quality 
factor  k and  then  express  the  two  corrections  as 


= kb  SSt  (43) 

and 

SC  ^ 

ACd  = k&  -f-  (44) 

V 

where  the  number  6 may  be  determined  from  the  literature  on  lift  inter 
ference  of  closed  tunnels. 

Solution  and  Interference  quality  factors  for  tunnels  with  cross 
sections  other  than  circular.-  The  solution  of  the  wall -interference 
problem  for  tunnels  having  cross  sections  other  than  circular  may  be 
obtained  by  following  the  method  of  paragraph  14.6,  reference  10.  This 
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method  requires  that  a function  z = f (£)  be  found  which  will  conformally 
transform  the  interior  of  the  tunnel  cross  section  in  the  z-plane  into 
the  interior  of  the  unit  circle  |£|  = 1 in  the  £-plane.  It  is  also 
necessary  that  f'(£)  does  not  vanish  or  became  infinite  in  the  unit 
circle  | £ | =1.  The  complex  velocity  of  the  tunnel  in  the  z-plane  is 
then 


u 


iv 


dw  dj, 
d f dz 


(45) 


In  the  function  dw/d£  of  equation  (45),  the  0n' s which  determine 
the  slot  edges  of  the  £ -plane  are  the  transformed  values  of  the  slot  edges 
in  the  z-plane  and  the  values  of  l in  the  £-plane  are  also  determined 
from  the  transformations  of  the  points  b in  the  z-plane.  Once  these 
values  are  determined  for  the  transformed  tunnel  in  the  £ -plane,  the 
velocity  field  dw/d £ in  the  £ -plane  may  be  computed.  The  velocities 
in  the  z-plane  may  then  be  computed  by  equation  (45). 

The  interference  and  the  quality  factor  may  now  be  deduced  from 
equation  (45).  Subtracting  the  free  field  frcm  the  complex  velocity 
given  in  equation  (45)  gives  for  the  interference  velocity: 


dwi  _ dw  d£  ibT 

dz  = d?  dz  + „(z2  . b2) 

The  interference  for  any  tunnel  can  be  determined  from  this  equation; 
however,  in  many  cases,  the  use  of  a quality  factor  may  be  more  con- 
venient. The  quality  factor  for  this  class  of  tunnels  may  be  written 


dw  dC  + ib  r 

k . dUz  "42-*2)  (lt7) 

dwc  dj;  ibr 

d"r  dz  + z2  _ b2) 

where  dwc/d£  represents  the  complex  velocity  of  the  transformation  of 
the  closed  tunnel  in  the  ^-plane. 

A simplification  of  equation  (47)  may  be  made  in  case  the  transforma- 
tion f(£)  may  be  approximated  by  ct;  for  points  near  the  center  of  the 
tunnel.  If  b is  sufficiently  small,  its  value  in  the  £-plane  may  there- 
fore be  represented  by  b = cl,  where  l represents  the  point  which 
locates  the  transformation  of  the  vortex  in  the  (>plane.  Substituting 
these  approximations  into  equation  (47)  gives  for  the  quality  factor: 
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dv  l icir 

k - ^ c rt( c2£2  - c2!2) 

dwc  £ + ictr 

d?  * „(c2t2  . c2j2) 

Equation  (48)  shows  that  the  quality  factor  near  the  center  of  the 
original  tunnel  in  the  z-plane  is  approximately  equal  to  the  quality 
factor  of  the  corresponding  circular  tunnel  in  the  £ -plane.  Thus,  to 
obtain  the  interference,  only  the  quality  factor  of  the  circular  tunnel 
needs  to  be  computed,  provided  the  interference  of  the  original  closed 
tunnel  is  known  and  the  approximation  z - cl  is  valid  throughout  the 
region  in  which  the  model  is  located. 


dw 


nr 


_ at  „(t2  - i2) 
Sis.  * il  l 


(48) 


dr,  *(r2  . i2) 


Applications  of  Theory  to  Various  Tunnels 

Two  slots  located  symmetrically  across  the  x-axis  (fig.  2(a)).- 
This  class  of  tunnels  has  the  symmetry  of  slot  location  that  is  treated 
under  case  IV,  so  that  its  complex  velocity  field  may  be  expressed  by 
the  product  of  equation  (34)  and  expression  (35d)  where  s takes  the 
values  0 and  1.  Since  the  integrals  across  the  two  panels  are  identically 
zero,  the  matrix  of  the  integrals  has  no  meaning  and  hence  is  not  used. 
With  these  considerations  the  complex  velocity  may  now  be  written 


dw 

dz 


ibr(l  - b^/b1*  - 2b2cos  20!  + 1 B2(z^  - l) 
F:- = (49) 

jt(z2  - b2)(l  - z2b2)y(z^  - 2z2cos  20^  + l)B2(b^  - l) 


The  wall  interference  is  determined  by  subtracting  the  complex 
velocity  for  the  free  field  from  equation  (49) , or 


dw^ 

dz 


ibr 

(1  - ,4), 

/b^  - 2b2cos  20^  + 1 

«U2  - b2) 

(l  - z2b2! 

)j - 2z2cos  20^  + 1 

(50) 


The  quality  factor  k is  determined  by  dividing  equation  (50)  by  the 
closed-tunnel  interference  which  is 


ibr  (l  - b2)(l  + z2) 
z2  - b2  _ 1 - z2b2 


(51) 
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or 


(l  - z^){\ 


b^  - 2bccos  20-^  + 1 


k = 


(l  - z2b^)  |/z4  - 2z2cos  20]_  + 1 

(l_-  ^X1  + z2)  . i 

(l  - z2b2) 


- 1 


(52) 


This  function  may  be  written 


k = 


(l  - z4)^b4  - 2b2cos  20-j_  + 1 - (l  - z2b2)^z4  - 2z2cos  2 6^  + 1 
(z2  - b2)  ^ z 4 - 2z2cos  20 + 1 


(53) 


The  effects  of  the  slots  on  the  interference  along  the  span  of  the  model 
are  obtained  by  substituting  x for  z;  then. 


k = 


(l  - x4)/b4  - 2b2cos  20x  + 1 - (l  - x2b 2)/x4  - 2x2cos  Z01  + 1 ^ 

(x2  - b2)^x4  - 2x2cos  20-^  + 1 


If  x and  b are  sufficiently  small  so  that  the  approximation 


& 


2x2cos  29j_  + 1 = 1-  x2cos  20^  + — sin220^ 

2 


is  valid,  then 


2 cos  20t  - 2x2  - (x2  + b2)sin220n 
k = i i (55) 

2 - 2x2cos  20-j_  + x4sin220]_ 

The  quality  factor  is  seen  to  be  insensitive  to  values  of  either  x or 
b,  so  that  for  this  tunnel  with  relatively  small  models  (b  = 0.25  or 
less)  the  effect  of  slotting  the  tunnel  on  the  lift  of  any  model  may  be 
obtained  by  multiplying  the  closed-tunnel  interference  for  that  model  by 
the  slotted-tunnel  quality  factor. 
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Pistolesi  (ref.  7)  shows  that  the  interference  quality  factor  for 
this  tunnel  with  a vortex  doublet  in  the  center  is,  in  the  notation  of 
this  paper,  cos  20]_.  Equation  (55)  reduces  to  the  same  value  when  x 
and  b are  set  equal  to  zero. 

Two  slots  located  symmetrically  across  the  y-axis  (fig.  2(b)).- 
This  class  of  tunnels  has  the  symmetry  of  Blot  location  that  is  treated 
under  case  III  (see  eq.  (35c)).  Thus,  its  complex-velocity  field  may  be 
expressed  by  the  product  of  equation  (34)  and  expression  (35c),  so  that 
it  may  be  written 


where  Aq  and  A2  are  determined  from  the  matrix  and  are 


Aq  = 


f0l  cos  20  d0 

J-&1  (l  - 2b2cos  20  + b4)^cos20  - cos20q 


(57) 


(58) 


If  z and  b are  considered  small  the  complex  velocity  of  this 
field  may  be  approximated  by 

dwl  _ ibT 
dz  rt 

from  which  the  quality  factor  becomes  approximately 


2Aq)  + (l  - b4)  (A2  - 2Aqz2)cos  20-l 
(A2  - 2A0b2)(l  - z2b2) 


[A2b2  - 2Aq)  + (l  - b4) (A2  - 2ApZ2)  cos  2 0]_ 


A2  - 2Aob‘ 
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and 


A3  = 


cob  e d0 


2b2cos  20  + +)  n /cos20  - cos20 
g=l 


g 


(63b) 


If  the  tunnel  contains  eight  slots  with  panels  intersecting  both 

axes. 


-ibr(l  - b*)  - 2b2cos  20g  4-  1 

dw  g=l 

— n 

1 - b*)(i . a*)  n | 2z^cos  20g  + 1 


A]_z^i 

(«8  ♦ 4) 

- V2' 

[*6  ♦ 

1 + a5i 

(zl°  + l) 

A^1 

[b2  + l] 

• - A3b2! 

[b6  + l)  + a5I 

[b“  + l) 

(64) 


where  Ai,  ky  and  A5  are  equal  to  the  determinants  remaining  when 

the  first,  second,  and  third  columns,  respectively,  are  removed  from  the 
following  matrix: 


'0 


1 cos  0 d0 


J-0n 


t(e) 


r 

Je0 


3 cos  0 d0 


f{9> 


®1  cos  30  cL0 

-e±  *'<«> 


'03  cos  30  d0 

0o  tie) 


P®1  cos  50  d0 

lex  “*(*> 


^3  cos  50  d0 

02  f(9) 


(6  5a) 


where 


f(e)  = (l  - 


2b2cos  20  + b*)  n ^cos  20 
g=l 


CO8^0 


S 


(65b) 
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If  the  tunnel  contains  12  slots  with  panels  intersecting  both  axes, 

6 

-ibr(l  - b4)  IT  l*k  - 2b2cos  29g  + 1 

dw  _ g^l x 

it  (z2  - b2)  (l  - z2b2)  _[  J_  ^z^  - 2z2cos  2 0g  + 1 

g=l 


A-^z^(z2  + l)  - A^z^z^  + l)  + A^z2^1^  + l)  - + .l) 

Axb6(b2  + l)  - A3b^(b6  + l)  + A5b2(b10  + l)  - A7(blJ+  + l) 


where  Aj_,  A^,  A5,  and  Ay  are  equal  to  the  determinants  remaining 

when  the  first,  second,  third,  or  fourth  columns,  respectively,  are 
removed  from  the  following  matrix: 


cos  0 d0  01  cos  3 0 d0 

.ei  He)  J-9l  ~ He) 


cos  5 6 d8 

f(e) 


1 cos  70  d0 


U f(0) 


03  cos  0 d0  f 03  cos  30  d0  f03  cos  50  d0  f03  cos  70  d0 


f(0) 


f(e) 


(67a) 


P05  cos  0 d0  f05  cos  30  d0  05  cos  50  d0  05  cos  70  d0 

J04  f(0)  J04  f(0)  J04  f(0)  J0I4.  f(0) 


where 


f (0)  = (l  - 2b2cos  20  + b^)  fl  /cos20  - cos^ 

g=l 


(67t) 


If  the  slots  are  of  equal  width  and  evenly  spaced,  equations  (62),  (64), 
and  (66)  can  be  simplified  by  substituting  the  terms 


1 - 2zmcos  20i  + zl 
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where  Ao,  A2,  and  A4  are  evaluated  from  the  following  matrix  by 
the  same  methods  used  in  equation  (64): 


r«i 

d0 

P91  cos  20  d0 

~\ 

J-0 1 

m 

J-^  He) 

U~ 

®1  cos  40  d0 

" i(e  ) 


i+01  d9 


- a 
3 ' 0l 


t(e) 


3-°l 


cos  29  d0 

f(8) 


f-ex 


(69a) 


COG  49  d9 

f(e) 


where 


f (0)  = (l  - 2b2cos  20  + bi+)^cos230  - cos2301 
For  case  IV,  slots  on  the  x-axis, 


(69b) 


dw 

dz 


-ibr( 

1 - th)\J-b12  - 2b6cos  60x  + 1 

B222 

r^i 

r-H 

1 

CO 

N 

1 

rH 

1 

-4 

*(*2  - b2 

) ^1  - z2b2)^z^2  - 2z8cos  602.  + 1 

B2b2(b^  - l)  - B4(b8  - 

■if 

(70) 


where 


B4  = 


3 91 


sin  20  d0 


(l  - 2b2cos  20  + b^)^cos^30i  - 008^30 


(71a) 


and 


rf -0i 


b2  = 


sin  40  d0 


(l  - 2b2cos  20  + b^O^cos^B],  - cos230 


(71b) 
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The  interferences  are  determined  as  in  the  previous  sections. 

Equations  (68)  and  (70)  may  be  extended  to  10,  14,  . . . 4n  + 2 
slots  in  the  same  manner  that  equation  (62)  was  extended  to  include  8 
and  12  slots. 


Single  slot  located  symmetrically  across  the  x-axis.-  Since  this 
tunnel  contains  an  odd  number  of  slots,  summation  (17 ) rather  than  (ll) 
must  be  used  in  equation  (14),  which,  when  expressed  for  a single  slot, 
becomes 

1 

z XZ  ( an  cos  — ^ + bn  sin 

dw  n=0  ' 2 

— 

(Z2  . -b2)(l  - z2b2)/  cos  -Q  - COS  0-J. 

Equation  (72)  may  be  rewritten,  with  several  modifications  after  sub- 
stituting z = e*1®,  as 


2n  + 1 \ 

T-* 

— — — (72) 


dw 

dz 


aQZ(l  + z)  + a-^(l  + z2)  + i|bpz(l  - z)  + bi(l  - z2Tj 
fz  (z2  _ b2)(l  _ z2b2)^l  - 2z  cos  0^  + z2 


(73) 


The  constants  in  this  case  must  be  evaluated  from  the  circulation  about 
the  poles  at  b and  -b  and  the  nonexistence  of  sources  at  those  points, 
as  the  potential  condition  is  automatically  satisfied  because  the  poten- 
tial is  constant  across  a single  slot.  The  circulation  and  nonexistence 
of  sources  may  be  expressed  about  the  pole  b as 


T + ia  = 


2«i^ob(l  + b)  + ai(l  + b2)  + i[bob(l  - b)  + bi(l  - b2T|| 

\jz  2b (l  - b^l  - 2b  cos  0X  + b2 


(74a) 


and  about  the  pole  -b  as 

2ni^ao(-b)(l  - b)  -*•  ax(l  + b2)  + i[b0(-b)(l  + b)  + bx(l  - b^)]} 


T + ia  = 


^2(-2b)(l  - b^}  ^1  + 2b  cos  + b2 


(74b) 


If  there  are  to  be  no  sources  at  the  poles,  a must  equal  zero.  It 
may  be  seen  from  inspection  that  ag  and  a^  are  equal  to  zero  and  that 
bo  and  bi  have  solutions  other  than  zero.  When  these  solutions  for 
bo  and  bi  are  substituted  into  equation  (73),  it  becomes 


CONFIDENTIAL 


L 


NACA  FM  L53A26 


CONFIDENTIAL 


33 


+■*  W 

II 


CvP  $ 

<<?  * 
M es 


CONFIDENTIAL 


This  factor  may  be  checked  hy  letting  01  = 0°,  which  represents  a closed  tunnel,  for  which 
may  be  shown  to  be  equal  to  +lj  and  if  0i  = l80°,  which  represents  the  open  tunnel,  then  ] 
be  shown  to  be  equal  to  -1. 
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Equation  (77)  may  be  considerably  reduced  by  using  the  following 
approximations,  in  which  z is  assumed  to  be  small  with  respect  to  1: 


2z  cos  0^  + 


9 sin20i 

1 - z cos  0]_  + zc  — s 


(78) 


/l  + 2z  cos  02_  + z2  = 1 + z cos  02_ 


>in20i 


+ z 


(7y) 


with  similar  equations  for  the  radicals  containing  b.  After  multiplica- 
tion and  collection  of  terms,  equation  (77)  reduces  with  the  above 
approximations  to 


k = 


sin2©. 

/ b sin20-i  \ 

COS  0-, 

1 2 

z 

1 + 4 - 2 ] 

? sin20-i 

1 - z cos  02_  + — - — - 


(80) 


If  terms  of  the  order  z2  and  higher  are  eliminated,  equation  (80)  may 
be  approximated  with 


Equation  (8l)  shows  that  the  interference  of  this  tunnel  has  an  odd 
function  component  along  the  x-axis  which  will  produce  a variation  in 
the  effective  angle  of  attack  that  will  cause  the  model  to  have  a 
rolling  moment.  Since  extension  of  flow  fields  of  this  type  (an  odd 
number  of  slots  symmetrically  located  with  respect  to  the  x-axis)  to 
greater  numbers  of  slots  will  not  eliminate  the  odd  powers  of  z,  the 
rolling  moments  will  continue  to  exist.  It  is  therefore  to  be  concluded 
that  tunnels  containing  an  odd  number  of  slots  symmetrically  located 
with  respect  to  the  x-axis  introduce  an  extraneous  moment  into  the  data. 

Single  slot  located  symmetrically  across  the  y-axis.-  Since  this 
tunnel  also  contains  an  odd  number  of  slots,  summation  (17)  must  be  used 
to  insure  the  proper  slot  symmetry.  With  this  choice  of  functions, 
equation  (15)  may  be  written, 


dw 


1q  cos  I + a2_  cos  ^ + bo  sin  ^ + b]_  sin  ~ 


dz  ‘ (z2  - b2)(l  - z2b2) 


^sin  $ - sin  01 
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Since  z = e*^,  the  complex  velocity  may  be  rewritten  after  sub 
stitution  and  combining  various  terms  as 


dw  (l  + i)|aQz(l  + z)  + &i{l  + z3)  + bQiz(l  - z)  + b]_i(l  - z3) 


dz 


2(z2  - b2)(l  - z2b2)  |/l  - 2iz  sin  0^_  - z ^ 


(83) 


The  constants  aQ,  aj_,  bo,  and  b^  must  now  be  evaluated  from 

the  circulation  about  each  of  the  poles  and  the  condition  of  continuity. 
The  two  conditions  about  each  pole  may  be  expressed  as 


r + io  = 


2ni(l  + i) 


aQb(l  + b)  + ax(l  + b3)  + b0ib(l  - b)  + b^l  - b3) 

r ■ . ——————  ■ — - ■ .....  — — 

4b  (l  - - 2ib  sin  0^  - b2 


(84a) 


about  the  pole  b and 


2rti(l  + i)|-aob(l  - b)  + ax(l  - b3)  - b0ib(l  + b)  + bji (l  + b3)l 
- r + ia  — ■— — - — 

-4b(l  - b2)l/l  + 2ib  sin  0-.  - b2 

(84b) 

about  the  pole  -b.  The  strength  a must  equal  zero  in  order  to  satisfy 
the  condition  of  nonexistence  of  sources.  In  order  to  simplify  the 
rationalization  of  equation  (84),  the  following  substitutions  are  made: 


and 

Also  let 


and 


C + iD  = /l  - b2  + 2ib  sin  0X 
C - iD  = {l  - b2  - 2ib  sin  0]_ 

Ai  = &Qb(l  + b)  + a]_(l  + b3) 

A g = -a<)b(l  - b)  + ax(l  - b3) 

B1  = bQb(l  - b)  + bx(l  - b3) 

B2  = -b0b(l  + b)  + bi(l  + b3) 

2bf(l  - b^)^l  - 2b2cos  291  + b^ 


> 


(85) 
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The  four  simultaneous  equations  (84)  may  now  be  written,  after  ration- 
alizing the  denominator, 


= R.P . (l  - 
= I.P.(1  - 
= R.P. (l  - 
= I.P.(1  - 


i)(C  + iD)(Ax  + iBL) 
i)(C  + iD)(Ax  + iBx) 
i)(C  - iD)(A2  + iB2) 
i)(C  - iD)(A2  + iB2) 


which  may  be  restated  as 


= AX(C  + D)  + B-^C  - D) 

= -AL(C  - D)  + B!(C  + D) 
= A2(C  - D)  + B2(C  + D) 
= -A2(c  + D)  + B2(c  - D) 


If  A^  and  B^  are  determined  from  the  first  two  and  A2  and  B2 
from  the'  last  two,  it  may  be  seen  that 


"4  - A1  - “0^(1  +t>)  * a^l  + to3) 

2(CZ  ♦ D ‘) 

^ ~ - a2  - -aoE(l  - E>  + ail1  - t3) 


?A4s-:  &T  - ei  - v<i  - f) - ^l1  - *3) 


21  Cc  + D 


y = ®2  = + *)  + bi(i  + b3) 

2(C2  + D2) 
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These  equations  may  be  solved  for  a q,  a^,  bQ,  and  b^  giving 


a0  “ 


al  ~ 


A 

(c  + 

D)(l  - b3) 

- (c  - 

D)(l  + 

b3) 

C2  + 

D2 

4b  (l  - 

. bM 

A 

(c  + 

D)b(l  - b) 

+ (c  - 

D)b(l 

+ b) 

C2  + 

D2 

H 

1 

.bM 

A 

Ic 

+ D)(l  - b3)  ~ (C 

- D)(l 

+ b3) 

+ I 

4b(  1 

- b4) 

A 

(C  + 

D)b(l  - b) 

+ (c  - 

D)b(l 

+ b) 

C2  + D2  4b(l  - b*0 


> 


(89) 


The  set  of  equations  (89)  may  be  simplified  to 


ao  = -b0  = 


2AI 


(d  - Cb3) 


4(c2  + D2)b(l  - b1*) 


= 2Ab(C  - Db) 

1 1 4(c2  + D2)b(l  - b^) 


a,  = bn  = 


r (90) 


Letting  hq  = -bo  and  a^  = b]_,  the  complex  velocity  equation  (83)  may 
be  rewritten 


dy 

2 

:»(« 

Iq  + alZ2) 

+ i 1 

(sqZ2  + aj 

dz 

2! 

(z2- 

bz; 

1(1  - .v) 

fir 

- 2iz  sin  0^  - z2 

(91) 
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Substituting  for  Eq  and  gives 


dv  _ 


2A  - Cb3  + b(C  - Db)z3  + i[b(C  - Db)  + (d  - Cb3)z^ 
4(c2  + D2)b(l  - b^)(z2  - b2)(l  - z2b2))/l  - 2iz  sin  0-^  - z2 


Nov  it  may  be  shown  that 


so  that 


dv  = _ £ 
dz  « 


where 


A = - ££l(l  - b^)(c2  + D2) 


,[(d  - Cb3)  + b(C  - Dbjz2!  + i jb(C  - Db)  + (d  - Cb3)z2| 

(z2  - b2)(l  - z^/T  - 2iz  sin  0i  - z2 


C + ID  = yl  - b*  + 2ib  sin  0X 


The  interference  velocity  is  obtained  by  subtracting  -ibr/jt(z2  - b2), 
so  that 

^ fb(C  - Db)  + (D  - Cb3)z2  -1 

ibr  >[d  - 0,3  + t(c  - Dblz2]  J ^ (gs) 

dz  n(z2  - b2)  b(l  - z2b2)^l  - 2iz  sin  0X  - z2 

Nov  divide  by  the  closed-tunnel  interference  to  obtain  the  quality  factor, 
so  that 

|b(C  - Db)  + (D  - Cb3) z2  - iz [d  - Cb3  + b(C  - Db)z2]  " 

|b(l  - z2b2)  {l  - 2iz  sin  0^  - z2 


n(z2  - b2) 


b(l  - z2b2Wl  - 2iz  sin  ©i  - z‘ 


-ibr 


n(z2  - b2) 
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b(C  - Db)  + (d  - Cb3)z2  - iz(p  - Cb3  + b(C  - Db)z2]|  3 
b(l  - - 2iz  sin  0]_  - z2  j 

k = i- wL  • 

b(z2  - b2)yl  - 2iz  sin  0-]_  - z2 

The  interference  quality  factor  at  the  center  of  the  tunnel  may  be 
examined  by  letting  z = 0;  then 

-bz 

From  the  evaluation  of  the  constants,  equation  (94)  > 


C = - 2b2cos  20i  + b^  cos  \|r 

D = ^ 1 - 2b2cos  20-l  + b^  sin  i|r 


when 


, _ 2b  sin  0i 

* = i tan"1  i (: 

2 1 - b2 

When  b is  small,  equation  (98)  may  be  simplified  by  first  approxi- 
mating C and  D.  These  approximations  axe,  for  small  values  of  b. 


2 

C = 1 - ^r(cos  20X  + sin20]_) 


D = b sin  0i 


When  the  values  for  C and  D given  in  equations  (102)  and  (103)  are 
substituted  into  equation  (98),  it  becomes 


k * sin  + 


008^01 
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This  equation  also  checks  with  the  results  given  in  reference  7> 
with  the  usual  change  in  sign  to  conform  with  the  notation  of  this  paper. 


Lift  interference  of  rectangular  tunnels.-  The  function  which 
transforms  a rectangular  tunnel  into  a circular  tunnel  is  given  in 
paragraph  Ih-. 8,  reference  10,  as 


5 = 


sn  y dn 

cn  ^2. 
2 


Xz 

2 


(105) 


where  sn  — , 
2 , 

and  x/2  is 


and 


defined  hy 


are  Jacobian  elliptic  functions  of  z 


X 

2 


K 

a 


E. 

h 


(106) 


where  K and  iK'  are  the  quarter  periods  and  a and  h are  the 
breadth  and  height  of  the  tunnel.  (See  pars.  14.7  and  14.8  of  ref.  10 
for  further  information  concerning  these  functions.) 

If  a and  h are  given,  then  K,  K',  and  m,  the  squared  modiolus 
(see  ref.  10,  par.  14.8),  are  uniquely  determined.  Once  the  constants  K, 
K' , and  m are  determined,  the  slot  location  and  the  half -span  length  Z 
may  be  computed.  First,  consider  slots  which  are  located  on  the  top  and 
bottom  of  the  tunnel  (see  fig.  lb).  In  reference  10,  paragraph  14.8,  it 

is  shown  that  the  top  of  the  tunnel  may  be  expressed  by  z = x + — or 

2 

Xz  = Xx  + iK' . It  is  also  shown  that  the  transformation  (105)  may  also 
be  expressed 


j,2  _ 1 - cn  Xz 


(107) 


1 + cn  Xz 

so  that  the  equation  for  the  top  of  the  tunnel  becomes,  in  the  £ plane, 


t2  - 1 - cn(Xx  + IK' ) 

1 + cn(Xx  + iK' ) 

From  reference  10,  paragraph  14.8,  cn(Xx  + iK')  is  equal  to 
which  may  be  called  in.  Now,  on  the  circle, 


(108) 

-mV2ds  iXx, 
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t 2 = e2iS  = 1 ~ 

^ 1 + ip 

and  if  the  real  and  imaginary  portions  are  separated, 


cos  20 


1 - u2  _ m - (is  Xx)2 
1 + m + (ds  Xx)2 


(109) 


(110) 


or 


sin  20 


2U _ 2ds  Xx  1 

2 ./-  , . ds  )ix 

+ M 1 + 


(HI) 


The  0n's  of  the  circular  tunnel  in  the  £ plane  are  determined  from 
the  above  equations  by  the  locations  of  the  slot  edges  xn  in  the  z-plane. 

If  the  slots  are  on  the  side  of  the  tunnel,  a similar  analysis 
starting  with  z=^+iy  or  Xz  = K + i Xy  will  show  that 


sin  20n 


2\fm  sd  iXy 
1 + m(sd  iXy)2 


(112) 


As  before,  equation  (112)  may  be  used  to  determine  the  0n's  when  the 
slot  edges  are  functions  of  y. 

The  value  of  l is  given  in  reference  10  as 


1 - cn  Xb 

l = 

1 + cn  Xb 


(113) 


Knowing  l and  the  0n's,  the  complex  velocity  of  the  tunnel  in  the 
£ plane  may  now  be  computed,  and  the  interference  determined  from 
equation  (16). 

The  calculation  of  the  interference  for  rectangular  tunnels  may  be 
simplified  by  applying  the  conclusion  following  equation  (18),  which 
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states  that  the  quality  factor  of  the  transformed  tunnel  is  approximately 
equal  to  the  quality  factor  of  the  original  tunnel,  provided  the  trans- 
formation function  can  be  approximated  with  z = c£.  Equation  (105)  is 
shown  in  reference  10,  paragraph  14.8,  to  be  approximately  equal  to 
£ = Xz/2;  hence  the  conclusion  following  equation  (48)  is  valid  for 
rectangular  tunnels.  Thus,  it  is  necessary  only  to  determine  the  quality 
factor  for  the  transformed  tunnel  and  applying  it  to  the  correction  for 
the  fully  closed  rectangular  tunnel  (see  ref.  10,  par.  14.8,  for  this 
correction)  to  obtain  the  correction  for  the  corresponding  slotted  rectan- 
gular tunnel. 

Circular  tunnels  having  general  slot  configurations.-  Equation  (l6) 
is  used  to  obtain  the  complex  velocity  for  a configuration  which  contains 
a general  slot  distribution.  Equation  (l6)  may  be  written 


dw 

dz 


m 

pz^  cos  i3  + bp  sin  -a)  (an  + i3n)cotn  ^ 

n=0 


(114) 


where  p has  the  same  definition  as  given  for  equation  (25).  If  z is 
substituted  for  i3  by  using  the  relation  z = eil3,  equation  (ll4)  becomes 


dw 

dz 


pz  (ai  Hr + lbl  Hr)  £ (°» + i(J»)  (*  rri)n 

ftm" 


(z2  - b2)(l  - z2b2)  nji  5-2-i  - cot  ^ 


(115) 


Equation  (115)  may  be  reduced  to 


— — I m 

p a-j_(l  + z2)  + ib^(l  - z2)  (on  + ipn)(z  - l)m“n(z  + l)nin 

dw  _ -*  n=0 

dz  2m 

i - M n 

6=1 


i(z  + 1)  - (z  - l)cot 


06 

T 


(116) 
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The  constants  a]_;  b]_;  a^,  a]_,  . . . a^;  and.  Pq>  3q»  . . . 3n 

are  determined  from  the  solution  of  the  set  of  simultaneous  equations  (117)* 
Once  the  constants  are  determined,  the  complex  velocity  may  be  determined, 
and  the  interference  may  then  be  computed  in  the  same  manner  as  was  used 
for  the  various  symmetrical  cases. 


RESULTS  AND  DISCUSSION 


Interference  quality  factors  for  several  circular  tunnels  with 
symmetrically  located,  evenly  spaced  slots.-  The  quality  factors  for 
circular  tunnels  containing  various  symmetrically  located,  evenly  spaced 
slots  of  equal  length  and  wings  of  very  small  span  (b  = 0)  are  given  in 
figure  3*  'The  curves  for  tunnels  that  have  panels  across  the  x-axis  and 
contain  4,  8,  or  12  slots  are  calculated  from  the  formula 


k 


cos ^5.  + jj  Q dQ 

l/cos^  — 0 - cos2 
V 2 2_ 

cos(|  - ly  9 d9 
j/cos^  — 0 - COS^ 

If  2 2 


(118) 


where  9 is  defined  in  figure  1 and  m is  the  number  of  slots.  Equa- 
tion (118)  may  be  derived  by  considering  the  value  of  equation  (34), 
expressed  for  symmetrically  located,  evenly  spaced  slots  of  equal  width, 
when  b and  z are  equal  to  zero.  It  is  also  the  negative  of  the 
relation  given  in  reference  7 for  the  same  tunnel  configurations,  the 
minus  sign  being  used  to  conform  with  the  notation  of  this  paper. 

An  analysis  of  figure  3 shows  that,  for  all  configurations  except 
the  one  with  two  slots  across  the  x-axis,  only  a small  percentage  of  the 
tunnel  wall  must  be  opened  in  order  to  obtain  no  interference  at  the 
center  of  a small  model,  and  that  the  amount  of  opening  required  rapidly 
becomes  smaller  when  the  number  of  slots  is  increased.  It  is  also  noted 
that  the  change  in  quality  factor  is  very  rapid  at  the  null-interference 
condition,  thus  making  an  accurate  estimate  of  the  interference  difficult. 

The  importance  of  slot  location  in  determining  the  quality  factor 
is  indicated  by  the  large  divergence  of  the  two  curves  for  tunnels  con- 
taining two  slots.  It  can  be  shown,  though,  by  comparing  the  quality- 
factor  functions  for  the  conditions  of  slots  across  the  x-  and  y-axes 
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and  panels  across  the  x-  and  y-axes,  that  as  the  number  of  slots  increases, 
the  quality  factors  will  approach  each  other.  Since  the  quality  factors 
approach  each  other,  it  may  be  expected  that  the  quality-factor  curves 
for  tunnels  containing  8 or  12  evenly  spaced  slots  of  equal  length  will 
be  approximately  correct  for  any  slot  location.  Thus,  the  quality  factor 
should  be  approximately  correct  even  though  the  model  is  rolled  in  the 
tunnel. 

Variation  of  quality  factor  with  span  of  model.-  The  effect  of  model 
span  on  the  quality  factor  k at  the  center  of  several  slotted  tunnels 
is  shown  in  figure  4.  Examination  of  the  curves  shows  that  the  quality 
factor  is  fairly  constant  over  an  appreciable  range  of  spans,  in  that 
it  does  not  vary  more  than  10  percent  for  spans  of  0.5  to  0.6  of  the 
tunnel  diameter  or  width.  The  quality  factor  for  the  tunnel  containing 
two  slots  is  seen  to  vary  more,  in  that  the  quality  factor  changes  0.10 
for  a change  in  the  span  from  zero  to  0.46  of  the  tunnel  diameter.  The 
quality  factor  of  the  single  slot,  however,  is  seen  to  remain  fairly 
constant  regardless  of  span. 

The  spanwise  variations  of  the  quality  factor  for  a model  with  a 
span  of  0.75  of  the  tunnel  diameter  is  shown  in  figure  5.  This  figure 
shows  that  the  quality  factor  for  the  tunnels  containing  8 or  12  slots 
is  about  0.3  larger  at  the  tips  than  it  is  at  the  center,  whereas  for 
the  tunnel  containing  two  slots  it  decreases  about  0.6.  These  changes 
indicate  that  the  spanwise  change  in  k may  depend  upon  the  boundary 
condition  at  the  intersection  of  the  wall  and  x-axis;  that  is,  if  a 
panel  intersects  the  x-axis,  k will  become  larger,  and  if  a slot  inter- 
sects the  x-axis,  k will  become  smaller.  Such  a tendency  would  be 
difficult  to  prove,  however,  so  that  an  analysis  of  the  variation  in  k 
for  any  particular  tunnel  would  require  a computation  of  k along  the 
span  for  that  model.  These  variations  in  k also  indicate  that  the 
lift  corrections  for  a model  with  a span  as  great  as  0.75  of  the  tunnel 
diameter  can  be  roughly  approximated  by  using  the  value  of  k at  the 
center  of  the  tunnel;  however,  a more  accurate  correction  would  involve 
the  use  of  an  average  value  for  k as  well  as  an  average  of  the  load 
over  the  span. 

A resume'  of  the  observations  made  from  figures  4 and  5 indicates 
that,  for  these  specific  slot  configurations  at  least,  the  quality 
factors  for  tunnels  containing  8 or  12  slots,  provided  they  are  symmet- 
rically located,  equally  spaced,  and  of  equal  length,  are  more  nearly 
constant  throughout  a greater  portion  of  the  central  region  of  the  tunnel 
than  are  the  quality  factors  in  tunnels  containing  two  slots.  The  fact 
that  the  quality  factor  is  moderately  constant  even  at  spans  of  0.75 
permits  making  the  correction  for  lift  interference  by  computing  the 
interference  of  the  closed  tunnel  having  the  same  cross  section,  and 
then  multiplying  that  interference  by  the  quality  factor  for  the  center 
of  the  slotted  tunnel. 
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Quality  factors  in  tunnels  containing  a single  slot.-  The  variation 
of  the  quality  factor  with  percent  opening  in  a tunnel  containing  a 
single  slot  symmetrically  located  with  respect  to  the  x-axis  is  shown 
in  equation  (80)  and  in  figure  6 to  be  somewhat  different  from  those 
previously  studied  in  that  k is  greater  on  the  half-span  of  the  wing 
which  points  toward  the  panel  and  is  smaller  on  the  side  which  points 
toward  the  slot.  This  variation  of  k means  that  a spanwise  variation 
in  angle  of  attack  exists  which  would  cause  the  model  to  roll.  Thus,  a 
tunnel  of  this  type,  that  is,  a single’  slot  symmetrically  located  with 
respect  to  the  x-axis,  has  interferences  which  are  more  difficult  to 
correct  for  because  of  the  introduction  of  an  unnecessary  rolling  moment 
into  the  measured  data. 

This  rolling  moment  does  not  exist,  however,  if  the  slot  is  symmet- 
rically located  with  respect  to  the  y-axis.  It  is  shown  in  equation  (97), 
though,  that  the  induced  velocities  will  not  be  normal  to  the  span  but 
will  have  a smal  1 component  along  the  span.  Such  a component  should  not 
affect  the  total  forces  seriously,  since  it  would  not  affect  the  total 
induced  velocity.  However,  the  effects  of  this  component  should  perhaps 
be  considered  in  the  treatment  of  load  distributions  along  the  span. 

Thus,  the  pecularities  of  the  interference  effects  of  tunnels  containing 
a single  slot  symmetrically  located  with  respect  to  either  axis  indicate 
that  its  effects  will  be  more  difficult  to  correct. 

If  the  equations  for  the  single-slot  case  are  extended  to  a larger 
odd  number  of  evenly  spaced,  equal -width  slots,  it  may  be  expected  that 
both  the  rolling  moment  and  the  cross  flows  will  become  smaller  because 
the  walls  create  a more  uniform  interference  field  at  the  model. 

Effects  of  compressibility  on  the  corrections.-  It  is  shown  in 
reference  11  that  the  lift  due  to  compressible  flow  can  be  corrected 
in  exactly  the  same  manner  as  though  the  flow  were  incompressible.  Thus, 
neither  the  k nor  the  6 of  equations  (43)  and  (44)  is  a function  of 
compressibility.  Since  the  arguments  used  in  reference  11  are  based  on 
subsonic  linearized  compressible  flow,  it  may  be  expected  that  correc- 
tions can  be  made  for  Mach  numbers  up  to  the  critical  or  slightly  higher. 

Applicability  of  the  theory.-  Several  of  the  differences  between 
the  idealized  and  the  actual  problem  are  those  due  to  viscosity,  which 
causes  a mixing  region  in  the  neighborhood  of  the  slot  rather  than  the 
assumed  constant  pressure  surface  which  divides  the  high-velocity  tunnel 
air  and  the  stagnant  tank  air.  Since  the  mixing  region  involves  various 
complicated  phenomena  such  as  turbulence,  velocity  gradients,  separation 
at  the  outer  slot  edges,  and  differences  in  boundary  conditions  in  dif- 
ferent slots  depending  upon  the  direction  of  flow  through  the  slots,  it 
is  very  likely  that  its  effects  on  the  lift  interference  will  have  to  be 
determined  from  analysis  of  experimental  data  (for  example,  ref.  12). 
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Another  effect  of  viscosity  which  becomes  important  if  the  slots 
are  narrow  and  deep  is  the  friction  on  the  air  as  it  flows  through  the 
slots.  Since  the  friction  reduces  the  amount  of  flow  through  the  slots, 
they  will  effectively  become  narrower,  so  that  the  quality  factor  may 
be  expected  to  become  larger. 

The  possibility  also  exists  that,  in  an  actual  tunnel,  the  pressures 
may  not  be  equal  in  all  the  slots  and  the  slot  pressures  may  not  be  equal 
to  the  pressure  in  the  tank.  It  would  be  expected  that  if  the  difference 
between  the  pressure  in  the  slot  and  that  in  the  tank  causes  more  flow 
through  the  slot  than  occurs  in  the  ideal  tunnel,  the  slot  will  effec- 
tively be  wider,  and  if  the  difference  decreases  the  flow,  the  slot  will 
effectively  be  narrower.  Thus,  the  actual  quality  factor  will  depend 
upon  the  effective  slot  widths  as  determined  by  the  pressure  increments. 

Other  differences  occur  because  a practical  tunnel  cannot  be  con- 
structed like  the  ideal  tunnel.  These  differences  involve  finite  slot 
lengths,  variable -width  slots,  and  lips  on  the  slot  edges.  In  con- 
sidering the  effects  of  finite  slot  lengths,  it  seems  reasonable  to  assume 
that  those  effects  should  be  no  more  serious  than  the  effects  of  a finite- 
length  open  tunnel.  In  reference  13,  it  is  shown  that  the  theoretical 
lift  corrections  for  an  infinitely  long  open  tunnel  are  adequate,  pro- 
vided the  model  is  located  a distance  of  at  least  one-half  the  tunnel 
height  from  the  entrance  and  exit  regions.  Therefore,  if  the  slot  con- 
figuration is  such  that  the  width  is  constant  over  a section  whose  length 
is  at  least  equal  to  tunnel  height,  the  theoretical  corrections  should 
be  adequate  even  though  the  slot  width  may  vary  considerably  outside  that 
region.  The  effect  of  lips  on  the  slots  may  be  shown  qualitatively  by 
comparing  the  pressure  gradients  of  the  two  types  of  slot  configurations, 
that  is,  with  and  without  lips.  Since  the  lips  confine  the  flow,  the 
pressure  gradient  in  that  configuration  is  less  steep  through  the  slot 
than  if  there  were  no  lips.  Therefore,  the  velocity  may  be  expected  to 
be  lower,  so  that  the  effective  width  of  the  slots  will  be  reduced, 
thereby  increasing  the  quality  factor. 


CONCLUDING  REMARKS 


An  analysis  of  the  equations  which  represent  the  interference  on 
the  trailing-vortex  system  of  a uniformly  loaded  wing,  due  to  wind- 
tunnel  walls  with  mixed  open  and  closed  boundaries,  has  shown  that: 

1.  Slot  openings  of  the  order  of  7 percent  for  four  evenly  spaced 
slots  of  equal  length,  and  less  for  larger  numbers  of  slots,  are  required 
to  reduce  the  interference  on  a lifting  model  to  zero.  The  zero- 
interference  quality  factor  is  critical  with  respect  to  the  percent  of 


CONFIDENTIAL 


48 


CONFIDENTIAL, 


NACA  RM  L53A26 


slot  opening  as  a small  change  in  slot  opening  will  cause  an  appreciable 
change  in  the  quality  factor  when  its  value  is  near  zero. 

2.  The  tunnels  which  contain  two  symmetrically  located  slots  showed 
quite  different  values  of  the  interference  for  different  slot  locations. 
The  differences  in  interference  became  smaller  as  the  number  of  slots 
increased. 

3.  In  the  tunnels  examined,  a region  existed  about  the  center  of 
the  tunnel  in  which  the  ratio  of  the  slotted-tunnel  interference  to  the 
closed- tunnel  interference  was  fairly  constant,  so  that  in  order  to 
obtain  the  corrections  for  the  effects  of  a slotted  tunnel  it  is  neces- 
sary only  to  multiply  the  closed-tunnel  interference  by  a constant. 

4.  The  region  in  which  the  ratio  of  the  slotted-tunnel  interference 
to  the  closed-tunnel  interference  is  reasonably  uniform  was  found  to  be 
larger  for  the  tunnels  containing  8 or  12  slots  than  for  those  containing 
2 slots. 


5.  An  examination  of  tunnels  containing  a single  slot  showed  that 
they  produced  a rolling  moment  or  a cross  flow  on  the  model.  These 
phenomena  may  be  expected  for  tunnels  containing  a larger  odd  number  of 
evenly  spaced,  equal-width  slots.  These  phenomena  should,  however, 
decrease  as  the  number  of  slots  increases.  These  phenomena  may  also  be 
expected  for  any  case  in  which  the  slots  are  asymmetrically  arranged 
with  respect  to  the  model  axes. 

6.  An  analysis  of  the  effects  of  compressibility  on  the  ratio  of  the 
interference  of  a slotted  tunnel  to  the  interference  of  a closed  tunnel 
shows  that  the  ratio  is  relatively  unaffected  by  compressibility  through- 
out the  subsonic  region. 
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Figure  3*-  Quality  factors  for  small-span  wings  (b  = 0)  in  circular 
tunnels  containing  various  numbers  of  evenly  spaced  slots  of 
equal  length. 
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Figure  4.-  Variation  of  the  quality  factor  at  the  center  of  the  model 

with  the  span  of  the  model. 
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Figure  5.-  Variation  of  the  quality  factor  along  the  span  in  several 
tunnels  containing  equally  spaced,  symmetrically  located  slots  of 
equal  length.  Span  of  model  Is  0.T5  of  tunnel  diameter. 
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Figure  6.-  Quality  factors  for  a tunnel  containing  a single  slot  and  a 
model  whose  span  equals  one-half  the  tunnel  diameter. 
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